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Partial Decode-Forward Relaying 
for the Gaussian Two-Hop Relay Network 

Jing Li and Young-Han Kim 


Abstract — The multicast capacity of the Gaussian two-hop 
relay network with one source, N relays, and L destinations 
is studied. It is shown that a careful modification of the partial 
decode-forward coding scheme, whereby the relays cooperate 
through degraded sets of message parts, achieves the cutset upper 
bound within (1/2) log N bits regardless of the channel gains 
and power constraints. This scheme improves upon a previous 
scheme by Chern and Ozgiir, which is also based on partial 
decode-forward yet has an unbounded gap from the cutset 
bound for L > 2 destinations. When specialized to independent 
codes among relays, the proposed scheme achieves within log N 
bits from the cutset bound. The computation of this relaxation 
involves evaluating mutual information across L(N + 1) cuts out 
of the total L2 N possible cuts, providing a very simple linear- 
complexity algorithm to approximate the single-source multicast 
capacity of the Gaussian two-hop relay network. 


I. Introduction 

Consider the Gaussian two-hop relay network with one 
source, N relays, and L destinations as depicted in Fig. [Q 
which can be viewed as a cascade of a broadcast channel 
(BC) from the source to the relays and multiple multiple 
access channels (MACs) from the relays to the destinations. 
The source node wishes to reliably communicate a common 
message to the L destination nodes with help of the N relays. 
The special case of L = 1, originally introduced by Schein and 
Gallager (T|, 0, is often referred to as the diamond network. 
The capacity is not known in general except for the trivial case 
of N = 1. 

The best known capacity upper bound is the cutset bound 
0, which is the maximum of the minimum mutual informa¬ 
tion across all possible cuts that separate the source and the 
destinations. There are several capacity lower bounds based on 
different coding schemes. The compress-forward scheme for 
the 3-node relay channel by Cover and El Gamal SI has been 
extended to relay networks in several forms, such as quantize- 
map-forward (QMF) by Avestimehr, Diggavi, and Tse 0, and 
noisy network coding (NNC) 0,0. The standard analysis 0 
shows that when specialized to our two-hop network model in 
Fig. □ these coding schemes achieve the cutset bound within 
O(N) bits for any channel parameters (recall that N is the 
number of relays). 
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Fig. 1. The Gaussian two-hop relay network. 


Recently, Chern and Ozgiir 0 provided a more refined 
analysis on the performance of NNC and showed that it 
achieves within (1/2) log 2 N(N + 1) bits from the cutset 
bound regardless of the number of destinations. In the same 
paper 0, Chern and Ozgiir extended the partial decode¬ 
forward (PDF) scheme for the relay channel by Cover and 
El Gamal 0 to the Gaussian diamond network (L = 1). In 
the PDF scheme by Chern and Ozgiir, the source broadcasts 
independent parts of the message to the relays, which in 
turn recover and forward their corresponding parts to the 
destination over the MAC. Thus, the Chern-Ozgiir scheme 
can achieve the rate characterized by the intersection of the 
BC capacity region and the MAC capacity region, which can 
be shown to be within log N bits from the cutset bound. When 
there are more than one destination node, however, the gap 
from the cutset bound becomes unbounded (8] Sec. VI]. 

In this paper, we develop an alternative extension of partial 
decode-forward that achieves the cutset bound within 4 log N 
bits for any number of destination nodes. In the proposed 
scheme, the relays decode for multiple message parts based 
on their respective decoding capabilities (as in the BC with 
degraded message sets 0) and forward these parts coopera¬ 
tively (as in the MAC with degraded message sets ED, ED- 
Thus, the proposed scheme achieves the rate characterized by 
the intersection of the capacity region of the BC with degraded 
message sets and the capacity regions of the group of multiple 
access channels with degraded message sets. 

Although this improvement may be viewed at first as an 
unnatural complication (except for the obvious benefit for 
achieving higher multicast rates with L > 2 destinations), it 
actually yields a simpler characterization of the achievable rate 
when independent Gaussian random codebooks are used at the 
relays, which yields a slightly looser but easier-to-compute 
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log N approximation of the capacity. A direct computation of 
the cutset bound as well as of the achievable rates for NNC 
and the Chern-Ozgiir PDF scheme requires evaluating mutual 
information across L2 N different cuts and then taking the 
minimum, which takes exponential time when directly com¬ 
puted. As an alternative to direct computation, approximate 
computation of the capacity (or the cutset bound) of the single¬ 
source single-destination relay network has been proposed by 
Parvaresh and Etkin m2 based on properties of submodular 
function minimization, which implies that the capacity of our 
two-hop network with L = 1 can be approximated within 2N 
in polynomial time of 0(LN 6 ) complexity (see also ITU ). 
In this paper, we refine and strengthen the Parvaresh-Etkin 
approximation result by showing that the achievable rate of our 
PDF scheme under independent codebooks involves evaluating 
only L(N + 1) cut rates. As a consequence, we develop an 
explicit algorithm to approximate the capacity as well as the 
cutset bound within log N with linear time complexity. 

Finally, we evaluate the performance of yet another variant 
of partial decode-forward for the two-hop relay network. 
Recently, Lim, Kim, and Kim developed distributed decode¬ 
forward, which generalizes partial decode-forward to general 
noisy networks for multicast M and broadcast CD. As in the 
case of noisy network coding, a naive analysis of distributed 
decode-forward results in an achievable rate within N/2 bits 
from the cutset bound. In this paper, we provide a refined 
analysis that establishes a gap of (log-ZV + 2) bits from the 
cutset bound. 

The rest of the paper is organized as follows. In the next 
section, we review basic facts on polymatroids. In Section UHl 
we formally define the capacity of the Gaussian two-hop relay 
network. In Section IIVI we review the cutset upper bound 
on the capacity, which will be benchmarked throughput. In 
Section[V] we review the Chern-Ozgiir partial decode-forward 
scheme for the Gaussian diamond network (L = 1). In Section 
IVII we present our coding scheme for the special case of the 
diamond network and then extend this result to the general 
/.-destination case. In Section IVIII we show the computation 
of the achievable rate of the relaxed version of our coding 
scheme involves linear complexity. In Section IVIIII and put 
forward the improved analysis of the performance of DDF. 
Finally, we conclude the paper. 

Throughout the paper, we mostly follow the notation in 
lf]~6l . In particular, we denote [1 : N ] := {1, 2, • • • , N}. The 
maximum of a finite set is denoted as J mayi := max(^7). A 
tuple of random variables is denoted as X(J) := (X ;j : j £ 
J). The Gaussian capacity function is defined as C(i) := 
(1/2) log(l + x). 


II. Mathematical Preliminaries 


Fet 0 : 2l 1:Ar l —> [0,oo) be a set function satisfying 

1) 0(0) = o, 

2) < 0(/C) if JCJC, and 

3) G/C) + U K)< f(J) + 0(/C). 


Then the polyhedron 

^(0) := |(ati,--- ,x N ) £ [0,oo)^ : 

<0G7), JC[l:iV] 

j£J 

is said to be a polymatroid (associated with 0); see, for 
example, D3- 

Example 1. For any random tuple (Xi,..., X,v. Y) such that 
Xi,... ,Xn are mutually independent, the set of rate tuples 
(Ri,... ,Rn) satisfying 

jej 

is a polymatroid EH Femma 3.1]. In particular, if Xj ~ 
N(0, Sj), j £ [1 : N ], and Y = iXj + Z, where 
Xi,..., X y and Z ~ N(0,1) are mutually independent, then 
the set of rate tuples (Ri ,... ,Rn) satisfying 


£ K iS°(l>) 

j£J \j£J J 

is a polymatroid. 

Example 2. Fet : [1 : N] —> [0,oo) be nondecreasing and 
define 0 : 2^ 1:N ^ —y [0, oo) by 


0G7) = 


0, 

$(Jmax), 


J = 0, 

otherwise. 


Then it can be readily shown that fX{(!>) is a polymatroid 
characterized by active inequalities 


k 

< 0([1 : A]) = *(*), ke[l:N}. 

3 =1 


In particular, for any random tuple (Xi, ..., Ay, Y), the set 
of rate tuples (Ri, ..., Rn ) satisfying 


k 

Y,Ri<i(x k -,Y |xf + i) 

i=i 


is a polymatroid. 

The following well-known result is pivotal in our discussion. 

Lemma 1 (Edmonds’s polymatroid intersection theorem |[T8l ). 
If J 2i (0) and <5^(0) are two polymatroids, then 


max I x i ' ( a ' 1 >''' j x n) 6 Z^(0) U ^(0) 
j=i 

= min [4>{J) + t/j(J c )]. 


III. Formal Definition of Capacity 

Recall the Gaussian two-hop relay network model depicted 
in Fig. [j] The received signals at the relays corresponding to 
the signal X from the source node are 

Yj = 9jX + Zj, j £ [1 :N], 
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where < 71 ,..., gjv are the channel gains from the source to 
relay nodes 1 through N, respectively, and Z\,...,Zm are 
independent N(0,1) noise components. We assume without 
loss of generality that 

I <?1 I > 1^21 > > | 9 n | • (1) 

Similarly, the received signals at the destinations correspond¬ 
ing to the signals Xi,..., X N transmitted from the relays are 

N 

Yd = g d jXj + Z d , d £ [1: L\, 

3 =i 

where g d j, j € [1 : TV], d € [1 : L], denote the channel 
gain from relay node j to destination node d, and Z \,..., Zl 
are independent N(0,1) noise components. The first (source- 
to-relays) hop of the network can be viewed as a Gaussian 
broadcast channel, while the second (relays-to-destinations) 
hop of the network can be viewed as multiple Gaussian 
multiple access channels. All nodes are subject to (expected) 
average power constraint P, and we denote by Sj = (jj P and 
Sdj = gljP the received signal-to-noise ratios (SNRs) at the 
relays and the receivers, respectively. 

We define a (2 nR ,n) code for a Gaussian two-hop relay 
network by 

• a message set [ 1 : 2 nR ], 

• an source encoder that assigns a codeword x N (m) to 
each message m £ [ 1 : 2 nR ], 

• a set of relay encoders, where encoder j £ [ 1 : N] assigns 
a symbol Xji{y l ~ l ) to each past received sequence y * -1 
for each transmission time is [ 1 : n ], and 

• a set of decoders, where decoder d S [1 : L] assigns 
an estimate rh d or an error message e to each received 
sequence y R . 

We assume that the message M is uniformly distributed over 
the message set. The average probability of error is defined as 
Pj"^ = P {Md 7 ^ M for some d S [1: L]}. A rate R is said to 
be achievable for the Gaussian two-hop relay network if there 
exists a sequence of ( 2 nR , n) codes such that lim n—>oo Pe ^ — 
0. The capacity C is defined as the supremum of all achievable 
rates. 

When N = 1, the capacity is 

C(S d i)j. 

For N > 2, however, no computable characterization of 
the capacity is known even when L = 1. In subsequent 
sections, we present bounds on the capacity and establish their 
closeness. 

IV. The Cutset Bound on the Capacity 

Since the network consists of two noninteracting channel 
layers, the cutset bound 0 on the capacity of a general noisy 
network can be simplified as 

C < Pcs 

:=supmin[l(X-Y(J c )) + I(X(J)-,Y d \X(J c ))\, (2) 

F 


C = min|C(S , : 


where the supremum is over all joint distributions F(x)F(x N ) 
satisfying E(X 2 ) < P and E(Xj) < P, j G [1 : N], 
the minimum is over all d £ [1 : L] and J C [1 : N\, 
and J c denotes [1 : N] \ J. By the maximum differential 
entropy lemma (see, for example, lfl6l Section 2.2]), the supre¬ 
mum in (0 is attained by Gaussian X and jointly Gaussian 
(Xi,..., X\r). By switching the order of the supremum (over 
Gaussian distributions) and the minimum, the cutset bound is 
further upper bounded as 


Pcs< sup mmsup[l(X;Y(J c )) + I(X(J)-Y d \X(J c ))] 

F{x N ) d X F{x) 


= sup nun 

F(x N ) d X . 


c (E 5 j 


I{X{J)-Y d \X{J c )) 


(3) 



Note that direct computation of the cutset bound in 0 for a 
fixed distribution or its relaxation in 0 involves evaluation of 
the minimum rate over the combination of 2 N choices of J 
and L choices of d, that is, the total L2 N cuts that separate 
the source and the destinations. 


V. The Chern-Ozgur Partial Decode-Forward 
Scheme for the Gaussian Diamond Network 


In the partial decode-forward scheme by Chern and 
Ozgiir Q3 (see also fl9ll ). which was developed mainly for 
the case N = 1, the source node divides the message M 
into N independent parts Mi,..., M ; \- (rate splitting), relay 
j recovers Mj and forwards it (decode-forward), and the 
destination node forms the estimates of Mi,..., Mjv and thus 
of M itself; see Fig. [2] This scheme is implemented over two 
hops in a block Markov fashion, and the achievable rate can 
be characterized as 


f N 

Ppdf = max rrP,:(Pi,. 
*G=i 


., Rn) € ^bc FI ^mac r ■ 


(5) 


Here ^bc is the capacity region of the standard A’-receiver 
Gaussian broadcast channel with SNRs Si,..., Sn, that is, 



Fig. 2. The Chem-Ozgiir partial decode-forward coding scheme for L = 1. 
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the set of rate tuples (i?i,..., Rn) such that 


Rj < C 


a j Sj 


ELi a kSk + 1 


j€[l:N], 


( 6 ) 


for some (cci,..., ajv) satisfying aj > 0, j € [1 : TV], and 
SjLi a j = 1> which, by the BC-MAC duality l20l , can be 
written as the set of rate pairs (Ri,..., Rjv) such that 


’ ^C[l:7V], (7) 

jeer S'eJ 7 


VI. The Proposed Partial Decode-Forward Scheme 
We propose a modified version of the Chern-Ozgur partial 
decode-forward scheme as depicted in Fig. [3] Here, the 
relays recover degraded sets of the message parts in the 
natural order—recall the assumption on the channel gains 
in CD — say, relay 1 recovers (Mi,..., Mjv), relay 2 recovers 
(M2, ..., Mjv), relay 3 recovers (M3,..., Mjv), and so on. 
The relays then cooperatively communicate these message 
parts to each destination as in the multiple access channel 
with degraded message sets ED, ED- 


for some (/3 1; ..., /3n) satisfying /3j > 0, j £ [1 : N], 
and J2 f=i Pj = 1- In 0- ^mac is the capacity region of 
the standard 7V-sender Gaussian multiple access channel with 
SNRs Sn, ..., S'ijv, i.e., the set of rate tuples (f?i,..., Rn) 
such that 

^TR^cAsA JC[l:N\. 
j&J 7 

Note that the region ^mac is a polymatroid (cf. Example [I]), 
but the region ftf-nc is not in general. Consequently, the 
maximum sum-rate of the intersection of the two regions, 
characterized by (0. is rather cumbersome to calculate. Chern 
and Ozgiir set /3j = 1 /N in 0 to obtain a polymatroidal inner 
bound on &bc characterized by 

^[1:7V], (8) 

jeJ v jeJ 7 

Now by 0 and Edmonds’s polymatroid intersection theorem 
with 


<KJ) = c EE). 

v j&j 7 


the corresponding (lower bound on the) achievable rate is 


Rpdf > min [4>(J) + 4>(J C )\ 

J-C[1:A] L J 


= min 
JC[1:N] 



+ c (e e) 


(9) 


By comparing this rate with the capacity upper bound in 0, 
we observe that the gaps for the two terms, both due to the 
lack of coherent cooperation, are bounded uniformly as 

c f£ s i)- c G?E s .i) s 5 losW ' (10) 

\j £ JC / V j G J c / 

c ((Ev / E)V c (ee) <iiogiv. (id 

V jeJ 7 7 


In conclusion, the gap between the achievable rate of the 
Chern-Ozgiir partial decode-forward scheme and the cutset 
bound is upper bounded as 


ApDF := RcS — Rpdf < log N , 
regardless of Sj and Sik, j , k £ [1: N\. 



Fig. 3. The proposed partial decode-forward scheme for L = 1. 


A. The Diamond Network 

For simplicity of exposition, we first consider the case 
L = 1. The achievable rate of the proposed scheme can be 
characterized as 

Rpdf = max /Rj : (f?i,..., Rn ) £ -^bc ^ ^mac 1 > 
S=i J 

where •^BC is the capacity region of the standard A’-receiver 
Gaussian broadcast channel (BC) with degraded message sets 
and is the capacity region of the A’-sender Gaussian 

multiple access channel (MAC) with degraded message sets. 
Since the broadcast channel is degraded in the order of 1 —> 
2 —>•••—} N, 3%bc = &bc as in 0. The capacity region of 
the multiple access channel with degraded message sets ED, 
m consists of all rate tuples (Ri,..., Rn) such that 
k 

J2Rj<I(X k ;Yi \X?+i), ke[l:N], (12) 
3=1 

for some F(x N ) such that E(X?) < P, j £ [1 : N], Again 
by the maximum differential entropy lemma, there is no loss 
of generality in setting (X±, ..., Xn) to be jointly Gaussian 
in (IT2l) . 

In order to obtain a lower bound on Rp DF , we follow the 
same approach 0 , mu as reviewed in the previous section 
and use the polymatroidal inner bound on &bc in 0. As for 
.AO we note the region in ( 1 1 21 is a polymatroid for a 
fixed F(x n ); cf. Example|2] Thus, by Edmonds’s polymatroid 
intersection theorem with 

<KJ) = I(X J — ;ni*SL, +1 ), ( 13 ) 

«^ = c (v£ s ->)' 

V 3EJ 7 
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the achievable rate of the proposed scheme is lower bounded Consequently, we can lower bound the achievable rate of the 
as scheme as 


Rpdf > sup min U)(J C ) + 4>(J)\ , (14) 

where the supremum is over all jointly Gaussian X N satisfy¬ 
ing E(Xy) < P, j G [1: N]. Since for each J C [1: TV] with 

\Jl nax — R 

+ <KJ) > m : k] c ) + <KJ) 

= f J ([k + l:N])+cj ) ([l:k}), 

the minimum in (ITTb is attained by J = 0 or J = [1 : A;] for 
some k. Thus, 


RpDF 


> sup min 
F fc£[0:iV] 


c 


N 

1 

N 

j=k +1 




+ I(X k -,Y 1 \X» +1 ) . 

(15) 


In comparison, by restricting J to be of the form [1 : k] in 
®, the cutset upper bound can be relaxed as 


Res < sup min 

F A:G[0:iV] 



+ i(x k ;Y 1 \x? +1 ) . 

(16) 


By comparing (fT5l > and (fl6l ). we establish the following. 


Proposition 1. The gap between the achievable rate of the 
proposed partial decode-forward scheme and the cutset bound 
is upper bounded as 

A'pdf := -Res ^ Rpdf A i log N, 
regardless of Sj and Sik, j , k 6 [1: N], 


B. The General Two-Hop Network 

The advantage of the modified partial decode-forward cod¬ 
ing scheme is fully realized when there are multiple destina¬ 
tions (L > 2), in which case the Chern-Ozgtir scheme has 
an unbounded gap from the capacity |8] Sec. VI]. Recall from 
Fig- [3]that in the proposed partial decode-forward scheme, the 
message parts are communication over a cascade of a BC (with 
degraded message sets) and multiple MACs with degraded 
message sets. The achievable rate can be thus characterized as 

Rpdf = max/ ^ Rj : (Ri,..., R_/v) € ^bc ^ ^mmac 1 > 

\7 = 1 J 

where M ac * s t ^ le set °f rate tuples (R \,..., Rn) such 
that 

k 

xk+1 ), kG [1:IV], (17) 

de[l:L\ 

3 = 1 

for some jointly Gaussian X N with E (Xf) < R, j G [1 : N], 
which is identical to the capacity region of the A’-sender in¬ 
state Gaussian compound MAC with degraded message sets. 

We can now proceed in the exactly same manner as in the 
single-destination case, except that in place of (fT3l> we have 
another polymatroid 

0(J) = 4 n G /( ^ maX;fd| ^- +l) - 


RpDF 


> sup min min 

F d£[l:L] fee[0:JV] 



+ I(X k ;Y d \X» +1 ) . 

(18) 


In comparison. 


Res 

< sup min min cf Sj ] + I(X k ] Y d \ X™, 

~ F de[l:i] fee[0:iV] V - / 

L '-j—k- )-l x 

This establishes the following. 


Theorem 1. The gap between the achievable rate of the 
proposed partial decode-forward scheme and the cutset bound 
is upper bounded as 

Apdf = Res — Rpdf < tj l°g N , 

regardless of the SNRs Sj and Sdk, j, k G [1: IV], d G [1 : L\, 
and the number of destinations L = 1,2,.... 


A few remarks are in order. 

1) When ^mac C ^bc (which is the case, for exam¬ 
ple, if |5 jv| > min d X]^Li I'M)’ the proposed coding 
scheme actually achieves the capacity 

In this case, the coding scheme simplifies to a simple 
decode-forward scheme, whereby every relay recovers 
the message M and coherently forwards it. 

2) At the other extreme, when ^bc — ^mmac (which is 
the case, for example, if |gi| < min,/ \g ( n |), the max¬ 
imum achievable rate of the proposed coding scheme 
is 

RpDF = C(Rl). 

Note that this rate is achieved trivially by using only the 
best relay (relay 1) and keeping the other relays idle, yet 
the gap from the capacity is no more than (1/2) log N. 
The performance difference between the Chern-Ozgtir PDF 
scheme and the proposed PDF scheme is best illustrated by 
the following example taken from J8] Sec. VI]. 

Example 3. Consider the Gaussian two-hop relay network 
with 2 relays and 2 destinations as depicted in Fig. [4] where 
the coefficients indicate the corresponding channel gains. The 
cutset bound is bounded as 


C(a 2 P) < Res < C((o + V^) 2 P), 

where the lower bound follows by setting X, X \, X 2 to be 
independent N(0, P) in (® and the upper bound follows by 
considering only the broadcast cut. The achievable rate of the 
PDF scheme by Chern and Ozgiir is 


Rpdf = C {aP), 
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These bounds yields a simple approximate expression for 
the capacity. 

Proposition 3. 


C = min min 
de[l:L] fce[0:AT] 


N 


S, 


j=k +1 


E* 

3=1 


dj 


±5 log N. 


Fig. 4. An example network. 


which has an arbitrarily large gap from the cutset bound as 
a —> oo. In comparison, the achievable rate of the proposed 
PDF scheme is lower bounded as 


R'pdf — C 


(a + a 2 )P 


which is within 1 bit from the cutset bound. 


VII. Linear-Complexity Capacity Approximation 

Computation of the achievable rate in (fl8l > requires max¬ 
imization over all Gaussian input distributions F. We now 
restrict the distribution to be independent and identically 
distributed Xj ~ N(0, P), j £ [1: N ]. This can be interpreted 
as a more practical coding scheme in which the relays use 
independent Gaussian codebooks and transmit codewords non- 
coherently. The achievable rate of the scheme is lower bounded 
by 


■ v j=k+1 7 v i=i 7 J 

(19) 

In comparison, starting with 0 and following the same 
argument as before, we can relax the cutset upper bound as 


> min min 
de[l:L]fce[0:IV] 


i?cs < min min 

de[l:L] fcG[0:iV] 


N 


s , 


j=k+l 



Thus, by ( TTTTb and (ITTI) . the capacity is approximated uniformly 
by log N. Moreover, the computation of ( fl9] > or (l20b involves 
computing Gaussian capacity functions for L(N +1) cuts, 
which is a significant savings from the directed computation 
of the cutset bound with all L2 N possible cuts as in 0. 

We summarize this result as follows. 


Proposition 2. The capacity of the Gaussian two-hop network 
is bounded as 

S,) + C[ZM ■ 

v i=fc+1 7 S=i 7 J 

S=fc+1 7 v j=i /J 

where the gap between the lower and upper bounds is no 
greater than log N for any Sj and S d k, j, k £ [1 : IV], d £ 
[1 : L\, and any L. Moreover, both bounds can be computed 
in 0(LN ) complexity. 


C > min min 

d£[l-.L] fce[0:A] 


C < min min 
de[l:L] k£[0:N] 


VIII. Distributed Decode-Forward 


In this section, we consider the distributed decode-forward 
(DDF) coding scheme in lfj~4l . which is an extension of partial 
decode-forward to general multicast networks. In particular, 
the rate achieved by DDF for our two-hop network is charac¬ 
terized M as 


I?ddf = sup min 

F d, 


I(X,X(J)-U(J c ),Y d \X(J c )) 


- Y, I(U k -,X,X N \Y k ) 

kGjc 


( 21 ) 


where the supremum is over all distributions of the form 
(UtiHik))F(x\x N )F(u N \x,x N ) satisfying E(X 2 ) < P 
and E(Xj) < P, j £ [1: N], By setting X and Xj to be i.i.d. 
N(0, P) and 


Uj =Yj — Zj + Z 3 , j £ [1: W], (22) 

where Zj ~ N(0,1), j £ [1 : IV], are independent of each 
other and of (X N , Y N ), it can be shown lfl4l that the gap 
between the achievable rate in ( OH and the cutset bound in 
0 is no greater than N/2. 

We now exploit the layered structure of the network to 
improve this O(N) gap to O(loglV). Following a similar (and 
in some sense dual) development for noisy network coding in 
|[8l , we set Zj ~ N(0, IV) in (02l >. Then, the first term of (OH 
becomes 


i(x,x(jy,u(j c ),Y d \x(j c )) 


( = } /(X; U(J c )) + I(X(J ); I X(J C )) 

= c (^E^-)+c(E£a 

v j&J c 7 y jeJ 7 


where (a) follows by the independence of (X, U N ) and X N 
and the layered structure of the network. For k £ [1: N\, each 
summand in the second term of (OH becomes 


I{Uk\X,X N \Y k ) = - log 


1 + (l + ~k)Sk 
1 + S k 


<o^g(l + - 


1 


< 


2 N' 


Hence, 


.Rddf > min 
d,J 


N 


j£J c 


E s d+c(E»«)- 




| J c 


2 N 


(23) 


Comparing this achievable rate in ( OH and the cutset bound 
in 0 establishes the following. 
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Proposition 4. The gap between the achievable rate of the 
distributed decode-forward scheme and the cutset bound is 
upper bounded as 

Addf = Res — -Rddf < log N + i, 

regardless of the SNRs Sj and Sdk, j , k £ [1 : N], d £ [1 : L], 
and the number of destinations L = 1,2,.... 


IX. Concluding Remarks 

Multiple coding schemes achieve the multicast capacity of 
the two-hop Gaussian network with one source, N relays, and 
L destinations within O (log IV), including: 

1) Noisy network coding (see fl9l Th. 3.1]) 

2) Distributed decode-forward (Prop. [4] in the current pa¬ 
per) 

3) Partial decode-forward (see jT9| Th. 3.3] for L = 1) 

4) Partial decode-forward with degraded message sets 
(Th. 1 in the current paper). 

Among these, the fourth scheme, which is the main contri¬ 
bution of the paper, achieves the tightest gap of (1/2) log N 
from the cutset bound. Moreover, a simple lower bound on its 
achievable rate can be expressed as the minimum of L(N 4-1) 
cut rates, providing a sharp approximation of the capacity 
that can be computed in O(LN) complexity. While it remains 
to be seen whether this linear-complexity approximation can 
be alternatively established via algebraic or combinatorial 
techniques, it is refreshing to note that the best computational 
result is obtained by a purely information-theoretic argument, 
based directly on a simple coding scheme. 
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